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Abstract. We exhibit a way to associate a quantum walk (QW) 
on the non-negative integers to any probabihty measure on the 
unit circle. This forces us to consider one step transitions that are 
not traditionally allowed. We illustrate this in the case of a very 
interesting measure, originally proposed by F. Riesz for a different 
purpose. 

For a review of Riesz's construction and its many uses, see [281 
[22l[TT]. For reviews of quantum walks, see [H IT2l [T3] . 



1. Introduction and contents of the paper 

The purpose of this note is to consider the probabihty measure 
constructed by F. Riesz, [18], back in 1918, and to study a quantum 
walk naturaUy associated to it. 

The measure on the unit circle that F. Riesz built is formally given 
by the expression 



d^{z) = n(i + cos(4^^))|^ = n(i + {^'^ + 



(1) 



2% 2Triz 

k=l k=l 



- j 

2'Kiz 

j=-oo 



E 



Here z = e^^ . If one truncates this infinite product the corresponding 
measure has a nice density. These approximations converge weakly to 
the Riesz measure. 
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The recent paper [B] gives a natural path to associate to a quan- 
tum walk on the non-negative integers a probability measure on the 
unit circle. This construction is also pushed to quantum walks on the 
integers. The traditional class of coined quantum walks considered in 
the literature allows for certain one step transitions and this leads to a 
restricted class of probability measures. 

In this paper we take the attitude that for an arbitrary probability 
measure a slightly more general recipe for these transitions gives rise 
to a quantum walk. The measure considered by Riesz falls outside 
of the more restricted class considered so far, and is used here as an 
interesting example. 

There is an obvious danger that having gone beyond the traditional 
class of quantum walks some of the appealing properties of these walks 
may no longer hold. From this perspective we consider the example 
discussed here as a laboratory situation where we will test some of these 
features. 

There is an extra reason for looking at this special example: Riesz's 
measure is one of the nicest examples of purely singular continuous 
measures. This means that the unitary operator governing the evolu- 
tion of the corresponding quantum walk has a pure singular continuous 
spectrum. Hence, Riesz's quantum walk becomes an ideal candidate to 
analyze the dynamical consequences of such a kind of elusive spectrum. 

We first show how to introduce a quantum walk given a probability 
measure on the unit circle and then we analyze in more detail the case 
of Riesz's measure, and give some exploratory results pertaining to the 
large time behaviour of the "site distribution" for this non-standard 
walk. 

We are grateful to Prof. Reinhard Werner for pointing out that 
F. Riesz actually started the infinite product ([T]) with k = 0. There are 
two well known references [221 E] that use the convention used here. 
Each choice has its own advantages as will be seen in section 5. 

For a review of Riesz's construction and its many uses, see [28] [22] 
[TT] . For reviews of quantum walks, see [U [121 [IS]- 

This paper will appear in the Proceedings of FoCAM 2011 held 
in Budapest, Hungary, to be published in the London Mathematical 
Society lecture Note Series. 

2. SZEGO POLYNOMIALS AND CMV MATRICES 

Let be a probability measure on the unit circle T = {z E C : 
\z\ = 1}, and L^^(T) the Hilbert space of /i-square-integrable functions 
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with inner product 

if,9) = / f{z)g{z)dp{z). 

For simphcity we assume that the support of fi contains an infinite 
number of points. 

A very natural operator to consider in our Hilbert space is given by 

(2) U,: Lj(T)^Lj(T) 

f{z) zf{z). 

Since the Laurent polynomials are dense in L^(T), a natural basis to 
obtain a matrix representation of f/^ is given by the Laurent polyno- 
mials {Xjjj^o obtained from the Gram-Schmidt orthonormalization of 
{l,z,z-i,zV-2,...}in Ll{T). 

The matrix C = {Xjy zXk)j^k=o ^f with respect to {xj} has the 
form 
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where pj = a/1 — and {aj}°^o ^ sequence of complex numbers 
such that \aj\ < 1. The coefficients aj are known as the Verblunsky 
(or Schur, or Szego, or reflection) parameters of the measure /i, and 
establish a bijection between the probability measures supported on 
an infinite set of the unit circle and sequences of points in the open 
unit disk. The unitary matrices of the form above are called CMV 
matrices, see [221 123| [27] . 

The problem of finding the sequence {aj} for a given measure p 
or, more generally, that of relating properties of the measure and the 
sequence is a central problem, see [22], where a few explicit examples 
are recorded. Even in cases when the measure is a very natural one, 
this can be a hard problem. Back in the 1980's one of us formulated 
a conjecture based on work on the limited angle problem in X-ray 
tomography. The same conjecture was also made in a slightly differ- 
ent context in work of Delsarte, Janssen and deVries. The conjecture 
amounts to showing that the Verblunsky parameters of a certain mea- 
sure are all positive. This was finally established in a real tour-de-force 
in [H]. 
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One of the results of this paper consists of finding these parameters 
in the case of F. Riesz's measure. In the process of finding these pa- 
rameters we will need to invoke some other sequences. Some of these 
will be subsequences of {aj}, and some other ones will only have an 
auxiliary role. We will propose an ansatz for the Verblunsky param- 
eters of the Riesz measure that have been checked so far for the first 
6000 non-null Verblunsky parameters. This is enough for computa- 
tional purposes concerning the related quantum walk. A proof of our 
ansatz deserves additional efforts. 

The decomposition of a measure dfi above into an absolutely contin- 
uous and a singular part can be further refined by splitting the singular 
part into point masses and a singular continuous part. The example of 
Riesz that we will consider later will consist only of this third type of 
measure, and it is (most likely) the first known example of a measure 
of this kind, built in terms of a formal Fourier series. For the case of 
the unit interval there is a construction of such a singular continuous 
measure in the classical book by F. Riesz and B. Sz-Nagy which is most 
likely due to Lebesgue. Notice that the method of "Riesz products" 
introduced in [IB] can be used to produce measures such as the one that 
lives in the Cantor middle-third set. However in this well as in 

the one due to Lebesgue, one loses the tight connection with Fourier 
analysis that makes the example of Riesz easier to handle. 

A very important role will be played by the Caratheodory function 
F of the orthogonality measure fi, defined by 

(4) F{z)= [l±^dfi{t), \z\<l. 

F is analytic on the open unit disc with McLaurin series 

oo „ 

(5) F{z) = 1 + 2 V/Z^.z^ /i, = / zU^x{z), 

3=1 ■^'^ 

whose coefficients provide the moments /ij of the measure /i. 

Another useful tool in the theory of orthogonal polynomials on the 
unit circle is the so called Schur function related to /i by means of F 
through the expression 

f{z) = z-\F{z)-l){F{z) + l)-\ \z\<l. 

Since the Schur function, and its Taylor coefficients will play such an 
important role, see [9], we will settle the issue of names by sticking to 
the name Verblunsky parameters for those that could be also called by 
the names of Schur or Szego. 
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These functions obtained here by starting from a probabihty mea- 
sure on T can be characterized as the analytic functions on the unit 
disk B = {z e C : \z\ < 1} such that F(0) = 1, Re F{z) > and 
1/(2;) I < 1 for 2; G D, respectively. 

Starting from /o = /, the Verblunsky parameters ak = /a:(0) can 
be obtained through the Schur algorithm that produces a sequence of 
functions {fk} by means of 



one can obtain a continued fraction expansion for f{z). This is called 
a "continued fraction-like" algorithm by Schur, [2T], and made into an 
actual one by H. Wall in [26j. See also [22]. We will illustrate the 
power of this way of computing these parameters by using it in our 
example to compute (with computer assistance, in exact arithmetic) 
enough of them so that we can formulate an ansatz as to the form of 
these parameters. 



We consider one-dimensional quantum walks with basic states \i) ® 
It) and \i) ® where i runs over the non-negative integers, and with 
a one step transition mechanism given by a unitary matrix U. This is 
usually done by considering a coin at each site i, as we will see below. 

One considers the following dynamics: a spin up can move to the 
right and remain up or move to the left and change orientation. A spin 
down can either go to the right and change orientation or go to the left 
and remain down. 

In other words, only the nearest neighbour transitions such that the 
final spin (up/down) agrees with the direction of motion (right/left) are 
allowed. This dynamics bears a resemblance to the effect of a magnetic 
interaction on quantum system with spin: the spin decides the direction 
of motion. This rule applies to values of the site variable i > 1 and 
needs to be properly modified at i = to get a unitary evolution. 





By using the reverse recursion 





3. Traditional quantum walks 
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Schematically, the allowed one step transitions are 



® It) 



\i) ® li) 



K + l)®lt) 
K - 1) ® li) 

K + l)®lt) 
K - 1) ® li) 



with amplitude 
with amplitude C21 

with amplitude 0*^2 
with amplitude C22 



where, in the case i = 0, the unitarity requirement forces the identifi- 
cation (g) II) = |0) (g) It). For each z = 0, 1,2, . . . , 
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is an arbitrary unitary matrix which we will call the i coin. 
If we choose to order the basic states of our system as follows 

(9) |0) ® It), |0) ® li), |1) ® It), |1) ® li), |2) ® It), |2) ® li), . 
then the transition matrix is given below 

\ 
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and we take this as the transition matrix for a traditional quantum 
walk on the non-negative integers with arbitrary (unitary) coins Ci as 
in g for z = 0, 1, 2, . . . . 

The reader will notice that the structure of this matrix is not too 
different from a CMV matrix for which the odd Verblunsky parameters 
vanish. This feature will guarantee that in the CMV matrix the central 
2x2 blocks would vanish identically. The CMV matrix should have 
real and positive entries in some of the 2x1 matrices that are adjacent 
to the central 2x2 blocks, and this is not generally true for the unitary 
matrix given above. In [U] one proves that this can be taken care of 
by an appropriate conjugation with a diagonal matrix. This associates 
a measure n on the unit circle with the above matrix U. Indeed, U 
becomes the matrix representation of the operator f/^, defined in ([2]), 
with respect to an orthonormal basis of Laurent polynomials Xj dif- 
fering only by constant phase factors e*^^ from the standard ones {xj} 
giving the CMV matrix. 
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In [0] one considers the case of a constant coin Cj for which the 
measure fi and the function F{z) are exphcitly found. In this case, 
after the conjugation alluded to above, the Verblunsky parameters are 
given by 

(10) a, 0, a, 0, a, 0, a, . . . 

for a value of a that depends on the coin, and the function F{z) is, up 
to a rotation of the variable z, given by the function 

z — z~^ — 2ilm a 

F z) -- 



y/{z - z-^y + A\a\^ - 2Re a 
The corresponding Schur function is the even function of z 



z^ 



2^2 



It is easy to see that, in general, the condition f(—z) = f[z) is 
equivalent to requiring that the odd Verblunsky parameters of should 
vanish. Traditional coined quantum walks are therefore those whose 
Schur function is an even function of z. In terms of the Caratheodory 
function the restriction to a traditional quantum walk amounts to 
F{-z)F{z) = 1. 



4. Quantum walks resulting from an arbitrary 
probability measure 

One of the main points of [6] was to show that the use of the measure 
/X allows one to associate with each state of our quantum walk a complex 
valued function in i^^(T) in such a way that the transition amplitude 
between any two sates in time n is given by an integral with respect 
to fi involving the corresponding functions and the quantity z". More 
explicitly we have 

(11) {^lU""]^) = [ z''^{z)iiz) dfx{z), 

where ipi^z) = J^ji^j-^ji^) the L'j^iT) function associated with state 
|\I') = '^j'ipjlj)- Here \j) is the j-th vector of the ordered basis con- 
sisting of basic vectors as given in (|9]), i.e. \j) stands for a site and a 
spin orientation, while Xj (z) are the orthonormal Laurent polynomials 
related to the transition matrix of the quantum walk. Similarly ipi^) 
is the function associated to the state 

This construction will now be extended to the case of any transi- 
tion mechanism that is cooked out of a CMV matrix as above. More 
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explicitly, we allow for the following dynamics 

with amplitude Pi+2Pi+3 
with amplitude Pi+iai+2 
with amphtude — Q;j+iaj+2 
with amphtude pj+2ai+3 

with amplitude — ai+2Pi+z 
with amplitude pj+ipj+2 
with amplitude — aj+ipi+2 
with amplitude — aj+2ai+3. 

The expressions for the amplitudes above are valid for any even i with 
the convention |— 1) ® \\) = |0) ® |t)- If i is odd then in every amplitude 
the index i needs to be replaced by i — 1. 

This generalization of the traditional coined quantum walks consists 
in adding the possibility of self-transitions for each site. One can, 
in principle, consider even more general transitions. As long as the 
evolution is governed by a unitary operator with a cyclic vector there 
is a CMV matrix lurking around. In our case the basis is given directly 
in terms of the basic states and there is no need to look for a new basis. 

It is clear that the main results in [B] extend to this more general 
case. If we have a way of computing the orthogonal Laurent polynomi- 
als we get an integral expression for the transition amplitude for going 
between any pair of basic states in any number of steps. 



N) ® It) 



\t) ® It) 



K + 1)® It) 
K-l)®li) 
K) ® It) 
IK)® It) 

^ + 1)8) It) 
^-1)® It) 

^) ® It) 
IK)® It) 



5. The Schur function for Riesz's measure 

From the expression for the Riesz measure given earlier we see that 
the expansion 

dz 



z'- 



2Tiiz 

j=-oo 

leads to the moments Hj of the measure p. Apart form the first one, 
Po = 1, if j 7^ can be written, in the necessarily unique form, as 

J = ±4^^! ± 4^=2 ± ■ ■ ■ ± 4*=^ ki> k2> ■■■ > kp>l, 

then 

p,- = 1/2^. 

For values of j that cannot be written in the form above we have pj = 0. 
In particular for j = 4^ we have fij = 1/2. 
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The moments of /i provide the Taylor expansion of the Caratheodory 
function 

oo 

and from this it is not hard to compute the first few terms of the Taylor 
expansion of the Schur function f{z) around z = 0. Indeed, from 

^12 ^20 ^44 ^48 ^52 ^60 

K J -r ^2 2 4 2 4 2 

we get that f{z) has the expansion 

•'^^^ ~ 2" ~ T ^ ~8~ ^ T6~ ~ 32 ~ ~64 128 256~ ^ " ' ' 

Only powers differing in multiples of 4 appear in the Taylor ex- 
pansion of both functions, F and /. This follows from the fact that 
diJ,{z) = dv{z'^) with v given by the same infinite product ([T]) as /i but 
starting at A; = 0, i.e. dfi^z) = {1 + {z + z~^)/2) dfi{z). From this 
we find that F{z) = Glz"^) and f{z) = z^g{z'^) where G and g are the 
Caratheodory and Schur functions of u respectively. 

It is now possible, in principle, to compute as many Verblunsky 
parameters for the function g as one wishes; they are given by the 
continued fraction algorithm given at the end of section |2} The first 
few ones are given below, arranged for convenience in groups of eight. 
We list separately the first four parameters. 
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-1/3 


5/8 


-1/13 


1/14 


-1/15 


-1/4 


-1/9 


1/10 


-1/11 


21/32 


-1/53 


1/54 


-1/55 


-3/52 


-1/49 


1/50 


-1/51 


5/56 


-1/61 


1/62 


-1/63 


-1/20 


-1/57 


1/58 


-1/59 


-11/48 


-1/37 


1/38 


-1/39 


-1/12 


-1/33 


1/34 


-1/35 


1/8 


-1/45 



It is clear that we get the Verblunsky parameters of / by introducing 
three zeros in between any two values above (a consequence of the 
argument z'^ in g above) and then shifting the resulting sequence by 
adding three extra zeros at the very beginning (a consequence of the 
factor z^ in front of g above), yielding finally the following sequence of 
Verblunsky parameters for /, where each row contains eight coefficients 
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starting with in the first row, ag in the second one, etc. 












1/2 











-1/3 











5/8 











-1/13 











1/14 











-1/15 











-1/4 
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21/32 
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-3/52 
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1/50 











-1/51 











5/56 











-1/61 



The non-zero Verblunsky parameters of / are given by 

= a4m-i, m = 1,2,3,... 

where the sequence {C,m} will be determined below. In fact it will be 
enough to determine the subsequence {^4+8n} since all the other values 
of C,m can be given by simple expressions in terms of these. 
The expression for these ^4+8n = «i5+32n is given by 

ai5+32n = -1/^n+l, 72 = 0,1,2,... 

for a sequence of integer values {An} to be described below. We will 
later give a different description of the complete sequence {^m} which 
makes clear what its limit points are and obviates the need to consider 
the subsequence {^4+8n}- 

We will first describe a procedure that allows us to generate the 
infinite sequence of integers {^nj^i of which the first ones are 

13, 53, 61, 37, 45, 213, 221, 197, 205, 245, 253, 229, 237, 149, 157, 133, 141, . 

Once this sequence is accounted for, i.e. if the Verblunsky parameters 
of the form ai5+32n are known, we will see that all the remaining ones 
are determined by simple explicit formulas in terms of {v4„}. For this 
reason we will refer to the sequence {An} to be constructed in the next 
section as the backbone of the sequence {a„} we are interested in. 

As we noted at the end of introduction, F. Riesz included the factor 
corresponding to = in the infinite product ([T]). That is, the measure 
considered by Riesz is the measure u giving our measure /i (starting the 
infinite product with /c = 1) by replacing z by z'^, and the corresponding 
Schur function is g. Therefore, the Verblunsky parameters {a^} that 
F. Riesz would have are those obtained deleting in the sequence {«„} 
the groups of three consecutive zeros, so = a4n+3 = C,n+i, and all of 
them should be computed from a§^_^ = — 1/A„. 
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The main difference between including the factor A; = in ([T| or 
leaving it out is the inclusion of many zeros in the list of Verblunsky 
parameters in the second case, which is the one we choose. This makes 
for a much sparser CMV matrix which is easier to analyze than it 
would be in the original case of F. Riesz. On the other hand his choice 
is better for computational purposes when, of necessity, one has to 
deal with truncated matrices. In Riesz's case there is more information 
packed in the same size finite matrix. This point is exploited in some 
of the graphs displayed at the end of the paper. 

6. Building the backbone 

Consider the sets Vj, j >0, defined as the ordered set of integers of 
the form 

((-2)^' - l)/3 + k2^+^ 

where k runs over the integers. As an illustration we give a few elements 
of the sets fo, fi, f2, ^"3, . . . , fio, namely 



Vo = .. 


.,-6,-4,-2,0,2,4,6,...; 


Vi = .. 


.,-13,-9,-5,-1,3,7,11,15,...; 


V2 = .. 


.,-23,-15,-7,1,9,17,25,...; 


V3 = ■■ 


.,-35,-19,-3,13,29,45,...; 


Vi = .. 


.,-91,-59,-27,5,37,69,101,...; 


V5 = ■■ 


. , -267, -203, -139, -75, -11, 53, 117, 181, . . . 


V6 = .. 


. , -491, -363, -235, -107, 21, 149, 277, 405, . . 


V7 = .. 


. , -1067, -811, -555, -299, -43, 213, 469, 725 


V8 = .. 


. , -1963, -1451, -939, -427, 85, 597, 1109, . . . 


Vg = .. 


. , -2219, -1195, -171, 853, 1877, 2901, . . . ; 


VlO = ■ 


. . , -5803, -3755, -1707, 341, 2389, 4437, .... 



These sets Vj , j > 0, are disjoint and their union gives all integers. 
A simple argument to prove this was kindly supplied by B. Poonen. 

We observe that dj, defined as the first positive element in the 
infinite set Vj (corresponding either to the choice = or A; = 1) is 
given as follows: if j = then do = 2 otherwise, for n > 1 we have 

d, = {{-2y - 1)/3 + (1 - (-iy)2^-. 

Define now, for n > 4, 

c„ = 8 + ((-2)"-^-l)2V3 
so that the values of C4, C5, Cg, C7, . . . are given by 

8, -24, 40, -88, 168, -344, 680, -1368, 2728, . . . 
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a sequence whose first differences are given by 

(-2r+\ n = 4,5,6,.... 

For each pair j, n, j > 0, n > 0, define wj^n as the number of positive 
elements in the sequence {vj} that are not larger than n. Notice that 
Yl'jLo '^3,n — ^ and that Wj^ is zero for all j. 

To be very explicit, we have for instance tyo,4o = 20, wi,4o = 10, 
W2,4o = 5, W3,4o = 2, ^^14 40 = 2, ^5,40 = 0, ^6,40 = 1, and all values of 
Wj^4Q after these ones vanish. 

It is possible to give an expression for Wj^n in terms of the sequence 
{dj} defined above. In fact one has 



n + 2^'+^ - dj 
2^+1 



where we use the notation \_x\ to indicate the integer part of the quan- 
tity X. 

Wc are finally ready to put all the pieces together and define, for 

n > 0, 

00 

(12) Un = ^Wj^nCA+j. 

3=0 

Notice that by definition this is a finite sum since, for a given n the 
expression „ vanishes if j is large enough. 

The reader will have no difficulty verifying that we get 

(13) A = 13 + 1*^-1, i>l. 
Recall that we have, for n > 0, 

1 



(14) Q;i5+32n = 

and, as mentioned earlier, we will see how all other values of aj can be 
determined from these ones. 



7. Building the sequence from its backbone 

We have observed already that the first non-zero Verblunsky pa- 
rameters of our measure are given by 

eta = 1/2, a'j — —1/3, an = 5/8, ctis = —1/13, .... 

We will see now that, for i > 15 and using the sequence {Ai\ built 
above, we have a way of computing the non-zero values of a^. Start 
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by observing that, after ctis = —1/Ai we get for values of j between 
j — 16 and j — 47 the following non-zero Verblunsky parameters: 

1 1 

"19 = :r"^ — ^, "23 



followed by 

_ 1 _ 1 _ 1 

""^^""A^' ""''-j;^' ""''-'A^' 

and finally 

The reader will have noticed that we did not give a prescription for 
a27 or for 0:43. This is done now: 

3 A2-A1 + 2 

We have seen that the non-zero values of aj for j in between j — 15 
and j = 47 are all obtained from the values of Ai and A2. We claim 
that exactly the same recipe apply for values of j in the range from 
16(2p — 1) — 1 and 16(2p -|- 1) — 1, with p > 1, namely we set 

1 1 

Q;i6(2p-1)-1 — Q;i6(2p+1)-1 — 

Zip 

and fill in the SEVEN non-zero values of aj in between these two by 
using expressions that are extensions of the ones above, namely, 

1 _ 1 

ai6(2p-l)+3 - Yj^' "l6(2p-l)+7 - ~2 + A ' 

_ 1 1 

"l6(2p-l)+15 — T) "l6(2p-l)+19 



1 1 



Q;i6(2p-l)+23 — 7:, Q;i6(2p-1)+31 



and, just as above, the missing Verblunsky parameters are given by 

«16(2p-l)+ll - -JTTY^ «16(2p-l)+27 - ^^^^ ^ _ 2 ■ 

8. A DIFFERENT EXPRESSION FOR THE VERBLUNSKY PARAMETERS 

The construction above gives as many non-zero Verblunsky param- 
eters for our measure as one wants, starting with aa, ay, an, aig, . . . 
for which we get the values 1/2, 1/3, 5/8, —1/13, .... 

In this section we give an explict formula for these Verblunsky pa- 
rameters in terms of a sequence of constants {Ki}, i — 0,1,2,..., 
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which are closely related to the sequence {Ai} introduced above. One 
of the advantages of this new expression is that the set of limit values of 
the sequence {aj} becomes obvious and is given by the union of three 
infinite sets, namely 

i = l,2,3,... , 

i = 0,l,2,3,... , 

z = 0,l,2,3,... . 

where the constants Ki are given by 
i^o = 3, 

X2,_i = 3A, X2i = 3(A-4), i = 1,2,3,.... 

One needs to add the limit points of the three infinite sets given above 
to get all limit points of the sequence {cxj}. 

We are now ready to give the alternative expressions for the non- 
zero Verblunsky parameters alluded to above, i.e. ^„ so that 

6 = 1/2, 6 = -1/3, 6 = 5/8,..., 

and in general a^rn-i = Cm- For this purpose we consider a disjoint 
union of the set of all non- negative integers into sets where the 
index n runs over the set 1,2,4,5,6,8,..., i.e. all positive n 7^ 3 
(mod 4). 

For each such n, define i?„ as the set of integers of the form 
1+4-^, p = 0,l,2,3,.... 

Once again, a simple proof of these properties of the sets S„ was sup- 
plied by B. Poonen. 

The sets break naturally into three classes, with n = (mod 4), 
n=l (mod 4) and n = 2 (mod 4) . We claim that 
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From the expressions above it follows that we have identified the 
limit values of the sequence {a^m-i}- The largest one is 2/3 and the 
lowest one —2/9. 

9. Some properties of the Riesz quantum walk 

Once we get our hands on the Verblunsky parameters correspond- 
ing to the Riesz measure we construct the corresponding CMV matrix 
and we can compute different quantities pertaining to the associated 
quantum walk. In the rest of the paper we choose to illustrate some of 
these results with a few plots. 

Figures [T] and |2] display the Verblunsky parameters themselves, key 
ingredients in the one-step transition amplitudes of the Riesz quan- 
tum walk. They show an apparent chaotic behaviour which is actu- 
ally driven by the rules previously described which generate the full 
sequence of Verblunsky parameters. This is in great contrast to the 
translation invariant case of a constant coin. 

Figures [3] and [5] display the Taylor coefficients of the Schur function 
for Riesz measure and the Hadamard quantum walk with constant coin 



on the non-negative integers. These coefficients have an important 
probabilistic meaning discussed in great detail in the forthcoming paper 
[9]: the n-th Taylor coefficient is the first time return amplitude in n 
steps to the state with spin up at site 0. 

The first time return amplitudes for the Riesz walk seem to fiuctuate 
in an apparent random way around a mean value which must decrease 
strongly enough to ensure that the sequence is square-summable, as 
Figure [6] makes evident. This is because the sum of the first return 
probabilities is the total return probability, which cannot be greater 
than one. Equivalently, any Schur function is Lebesgue integrable on 
the unit circle with norm bounded by one, and its norm is the sum of 
the squared moduli of the Taylor coefficients. 

In the Hadamard case the behaviour of the first time return am- 
plitudes is much more regular and the convergence to the total return 
probability, depicted in Figure |4| holds with a much higher speed. We 
should remark that the plot for the Hadamard walk picks up only the 
first 70 non-null coefficients, while the Riesz picture represents the first 
non-null 7000 coefficients. Hence, the differences between these two 
examples are not only in the more regular pattern that the Hadamard 
return probabilities exhibit, but also in the much higher tail for the 
Riesz return probabilities. 
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Figures [7] and |8] give the probability distribution of the random 
variable Xn/n, where X„ stands for the position (regardless of spin 
orientation) after n steps of the quantum walk started at position 
with a spin pointing up. The plots given here correspond to the value 
n = 800 for both, the Riesz quantum walk and the Hadamard constant 
coin on the non-negative integers. 

The figures show that, in contrast to classical random walks for 
which Xn behaves tipically as ^/n, the position in a quantum walk can 
grow linearly with n. Nevertheless, Figure |8] shows a striking behaviour 
of the Riesz walk compared to the more regular asymptotics of the 
Hadamard walk reflected in Figure [7j This should be viewed as a 
clear indication of the anomalous behaviour that can appear under 
the presence of a singular continuous spectrum. In particular, these 
results make evident that quantum walks with a singular continuous 
measure can not exhibit nice limit laws as other toy models do. For 
the case of translation invariant ones it is known that obey inverted 
bell asymptotic distributions (see for instance |il3j). 

These results should motivate a more detailed analysis of quantum 
walks associated with singular continuous measures. This could lead 
to the discovery of new interesting quantum phenomena. 
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Figure 1. Riesz's measure: The first 30000 non-zero 
Verblunsky parameters. 
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Figure 2. Riesz's measure: The non-zero Verblunsky 
parameters for indices between 30000 and 60000. 
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Figure 3. Hadmard's walk: The first 70 non-zero Tay- 
lor coefficients of the Schur function. 
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Figure 4. Hadamard's walk: The cumulative sums of 
the squares of the first 70 non-zero Taylor coefficients of 
the Schur function. 
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Figure 5. Riesz's measure: The first 7000 non-zero 
Taylor coefficients of the Schur function. 




Figure 6. Riesz's measure: The cumulative sums of the 
squares of the first 7000 non-zero Taylor coefficients of 
the Schur function in steps of 100. 
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Figure 7. Hadamard's walk on the non-negative inte- 
gers: The probabihty of X„/n for 800 iterations starting 
at |0) ® It). 
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Figure 8. Riesz's walk on the non-negative integers: 
The probabihty of Xn/n for 800 iterations starting at 
|0)®|t). 
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